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Abstract 

We define a stochastic model of a two-sided limit order book in terms of its key quantities best bid [ask] 
price and the standing buy [sell] volume density. For a simple scaling of the discreteness parameters, that 
keeps the expected volume rate over the considered price interval invariant, we prove a limit theorem. 
The limit theorem states that, given regularity conditions on the random order flow, the key quantities 
converge in probability to a tractable continuous limiting model. In the limit model the buy and sell 
volume densities are given as the unique solution to first-order linear hyperbolic PDFs, specified by the 
expected order flow parameters. 
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Introduction 

A significant part of financial transactions is nowadays carried out over an electronic limit order book 
(LOB) where unexecuted orders are stored and displayed while awaiting execution. From a mathematical 
perspective, LOBs are high-dimensional complex priority queueing systems. Incoming limit orders can be 
placed at many different price levels while incoming market orders are matched against standing limit orders 
according to a set of priority rules. Almost all exchanges give priority to orders submitted at more competitive 
prices (“price priority”) and displayed orders have priority over hidden orders at the same level (“display 
priority”). Orders with the same display status and submission price are usually served on a first-come-first- 
serve basis. The inherent complexity of limit order books renders their mathematical analysis challenging. In 
this paper, we propose a simple and transparent queueing theoretic LOB model whose dynamics converges 
to a coupled ODE:PDE system after suitable scaling. The limiting system can be solved in closed form. 

*Most of this work was done while the first author was visiting CEREMADE; grateful acknowledgment is made for hospitality 
and financial support. Both authors acknowledge financial support through the SFB 649 “Economic Risk” and Deutsche Bank 
AG. We thank seminar participants at various institutions for valuable comments and suggestions. 
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There is a significant economic and econometric literature on LOBs including Biais at al. [5], Easley and 
O’Hara [13], Foucault et al. Gloston and Milgrom (TB], Parlour [53], Rosu m and many others that 
puts a lot of emphasis on the realistic modeling of the working of the LOB, and on its interaction with 
traders’ order submission strategies. There are, however, only few papers that analyze order flows and the 
resulting LOB dynamics in a mathematically rigorous manner. Among the first was the one by Kruk [24] . He 
studied a queueing theoretic model of a transparent double auction in continuous time. The microstructure 
may be interpreted as that of a simple limit order book, if one considers a buyer as a buy limit order and 
a seller as a sell limit order. At the auction, there are TV S N possible prices of the security and thus 2N 
different classes of customers (N classes of buyers and N classes of sellers). Kruk established diffusion and 
fluid limits. The diffusion limit states that for N=2 possible prices, the scaled number of outstanding buy 
orders at the lower price and the scaled number of outstanding sell orders at the higher price converge weakly 
to a semimartingale reflected two-dimensional Brownian motion in the first quadrant. The fluid limit is such 
that various LOB quantities converge weakly to affine functions of time. 

Since information on the best bid and ask price and volume at different price levels is available to all market 
participants, it is natural to assume that the order arrival dynamics depends on the current state of the 
order book. The feature of conditional state-dependence was considered by Cont et al. [12], who proposed 
a continuous time stochastic model with a finite number of possible securities prices where events (buy/sell 
market order arrival, buy/sell limit order placement and cancelation) are modeled using independent Poisson 
processes. The arrival rates of limit orders depend on the distance to the best bid and ask price in a power- 
law fashion. The authors were able to show that the state of the order book, defined as a vector containing 
all volumes in the order book at different prices, is an ergodic Markov process. Using this fact, several 
key quantities such as the probabilities of a mid price move, a move in the bid price before a move in the 
ask price, or the probability of volume execution before a price move could be computed and benchmarked 
against real data without taking scaling limits. 

Cont and de Larrard m considered a scaling limit in the diffusion sense for a Markovian limit order market 
in which the state is represented by the best bid and ask price and the queue length, i.e. the number of 
orders at the best bid and ask price, respectively. With this reduction of the state space, under symmetry 
conditions on the spread and stationarity assumptions on the queue lengths, it was shown that the price 
converges to a Brownian motion with volatility specified by the model parameters in the diffusion limit. Very 
recently, Cont and de Larrard [10] studied the reduced state space under weaker conditions and proved a 
refined diffusion limit by showing that under heavy traffic conditions the bid and ask queue lengths are given 
by a two-dimensional Brownian motion in the first quadrant with reflection to the interior at the boundaries, 
similar to the diffusion limit result for N = 2 prices in Kruk [24] . 

In the framework analyzed by Abergel and Jedidi [1], the volumes of the order book at different distances to 
the best bid and ask were modeled as a finite dimensional continuous time Markov chain and the order flow 
as independent Poisson processes. Under the assumption that the width of the spread is constant in time, 
using Foster-Lyapunov stability criteria for the Markov chain, the authors proved ergodicity of the order 
book and a diffusion limit for the mid price. In the diffusion limit, the mid price is a Brownian motion with 
constant volatility given by the averaged price impact of the model events on the order book. 

In this paper, we prove a law of large numbers result for the whole book (prices and volumes). Specifically, 
we propose a continuous-time model of a two-sided state-dependent order book with random order flow and 
cancelation and countably many submission price levels. The buy and the sell side volumes are coupled 
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through the best bid and ask price dynamicslj We model the buy and sell side volumes as density functions 
in relative price coordinates, i.e. relative to their distance to the best bid/ask prices. Volumes at positive 
distances are defined as standing limit volume, observable in the order book (“visible book”), and volume at 
negative price distances corresponds to volume that would be placed in the spread if the next order would 
be a spread limit order (“shadow book”). The state of the model at any point in time is thus a quadruple 
comprising the best bid price, the best ask price, the relative buy volume density function and the relative 
sell volume density function. The state dynamics is then defined in terms of a recursive stochastic process 
taking values in a function space. 

In order to establish our scaling limit we use a mathematical framework similar to the recurrently defined 
semi-Markov processes in Anisimov ma and Gikhman and Skorokhod m, but for stochastic processes 
taking values in Banach spaces. When the analysis of the market is limited to prices as in e.g. Garman 
|16j , Bayraktar et al. [7] or Horst and Rothe |23j or to the joint dynamics of prices and aggregate volumes 
(e.g. at the top of the book) as in Gont and de Larrard [TUHH], then the limiting dynamics can naturally 
be described by ordinary differential equations or real-valued diffusion processes, depending on the choice of 
scaling. 

The analysis of the whole book including the distribution of standing volume across many price levels is much 
more complex. Osterrieder [28] modeled LOBs using measure-valued diffusions. Our approach is based on 
an averaging principle for Banach space-valued processes. The key is a uniform law of large numbers for 
Banach space-valued triangular martingale difference arrays. It allows us to show that the volume densities 
take values in and that the noise in the order book models vanishes in the limit with our choice of scaling. 

Our scaling limit requires two time scales: a fast time scale for cancelations and limit order placements 
outside the spread (events that do not lead to price changes), and a comparably slow time scale for market 
order arrivals and limit order placements in the spread (events that lead to price changes). The choice of 
time scales captures the fact that in real world markets significant proportions of orders are never executed. 
Mathematically, the different time scales imply that aggregate cancelations and limit order placements 
outside the spread between consecutive price changes can be approximated by their expected values. 

Our main result states that when limit buy [sell] orders are placed at random distances from the best ask 
[bid] price and the price tick tends to zero, order arrival rates tend infinity, and the impact of an individual 
placement/cancelation on the standing volume tends to zero, then the sequence of scaled order book models 
converges in probability uniformly over compact time intervals to a deterministic limit. The limiting model 
is such that the best bid and ask price dynamics can be described in terms of two coupled ODE:s, while 
the dynamics of the relative buy and sell volume density functions can be described in terms of two linear 
first-order hyperbolic PDE:s with variable coefficients. The PDE can be given in closed form. 

The remainder of this paper is organized as follows. In the first section we define a sequence of limit 
order book in terms of four scaling parameters: price tick, expected waiting time between two consecutive 
orders, volume placed/canceled and the proportion of order arrivals leading to price changes. In Section [2] 
we establish convergence of the bid/ask price dynamics to a 2-dimensional ODE; Section [3] is devoted to 
the analysis of the limiting volume dynamics. The uniform law of large numbers for triangular martingale 

^The coupling of the buy and sell sides through prices is essential for the limiting volume dynamics to follow a PDE. It is 
not essential for obtaining a scaling limit per se. Our mathematical framework is flexible enough to allow for dependencies of 
order flows on standing volumes, but that one would lead to a function-valued ODE as the scaling limit, rather than a PDE. 
The PDE-scaling is so much more transparent that it justifies, in our view, the restriction of order flow dependencies on prices. 

^We refer to Mandelbaum et al. m for general approximation results for queueing systems. 
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difference arrays as well as useful auxiliary results are proved in an appendix. 


1 A sequence of discrete order book models 


In electronic markets orders can be submitted for prices that are multiples of the price tick^ the smallest 
increment by which the price can move. In this section, we introduce a sequence of order book models for 
which we establish a scaling limit when the price tick and impact of a single order on the state of the book 
tend to zero, while the rate of order arrivals tends to infinity. 

The sequence of models is indexed by n G N. We assume that the set of price levels at which orders can 
be submitted in the n:th models is where Z denotes the one-dimensional integer lattice]! We put 

:= j ■ for j ^’L where is the tick size in the n:th model. 

The state of the book changes due to incoming order flow and cancelations of standing volume. The state 

(n) 

after fc G N such events will be described by a random variable ^ taking values in a suitable state space E. 
In the n:th model, the fc:th event occurs at a random point in time The time between two consecutive 
events will be tending to zero sufficiently fast as n —>■ oo. The state and time dynamics will be defined, 
respectively, as 

( 1 ) 
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random operators that will be introduced below. The conditional expected increment of the state sequence 
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1.1 The order book models 

In the sequel we specify the dynamics of our order book models. Throughout all random variables will be 
defined on a common probability space (H, P). 

1.1.1 The initial state 

The initial state of the book in the n:th model is given by a pair Aq^^) of best bid and ask prices 

together with the standing buy and sell limit order volumes at the various price levels. It will be convenient 
to identify the standing volumes with step functions 


(n) 


j=o " ' ■ “ " " 


j=0 


that specify the liquidity available for buying and selling relative to the best bid and ask price. The liquidity 
available for buying (sell side of the book) j G No ticks above the best ask price at the price level is 


'^j+i 


(n) 


(x)dx = ■ Ax^"i. 


®The assumption that there is no minimum price is made for analytical convenience and can easily be relaxed. 
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Figure 1: Snapshot of the book. 
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In order to conveniently model placements of limit orders into the spread, we extend q and q to the 
negative half-line. The collection of volumes standing at negative distances from the best bid/ask price is 
referred to as the shadow book. The shadow book will undergo the same dynamics as the standing volume 
(“visible book”). At any point in time it specifies the volumes that will be placed into the spread should 
such an event occur nexfl The role of the shadow will be further illustrated in Section 11.1.41 below; see 
Figures 3 and 4. 


Definition 1.1. In the n:th model the initial state of the book is given by a quadruple 

qM/ \ _ f tjM aM \ „(")/' 'lA 

*^0 — 1^0 I ^b,o Iv; ^s,0 i'lj 

where are the best bid/ask price and the step functions ® [0i oo) are to he 

interpreted as follows: 


-^ 6,0 ' 




f 

standing buy [sell] limit order volume density at price distance x 
below [above] the best bid [ask] price, for a; > 0 (visible book) 

< 

potential buy [sell] limit order volume density at price distance x 
above [below] the best bid [ask] price, for a; < 0 (shadow book). 


(3) 


Throughout, we shall use the notation / = 0{g) and / = o{g) to indicate that the function / grows 
asymptotically no faster than g, respectively that \f{x)/g{x)\ —>• 0 as a; —>• oo. With this, we are ready to 

"^Notice that the liquidity available for buying is captured by the sell side of the book and vice versa. 

®One has to specify the volumes placed into the spread somehow. Our choice of shadow books is one such way. The role of 
the shadow book will be further clarified in the following subsection when we define the impact of order arrivals on the state of 
the book. Its initial state is part of the model; future states will undergo the analogous dynamics as those of the visible book. 
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Figure 2: Left [right]: the buy [sell] volume density function in relative price coordinates. The colored area 
under the densities is the standing volume; the white is the shadow book. 


state our conditions on the initial states. In particular, we assume that the initial volume density functions 
vanish outside a compact price interval@ 

Assumption 1.2 (Convergence of initial states). The initial volume density functions vanish outside a 
compact interval [—M, M] for some M > 0 and and satisfies: 

± - u(”)(-)|oo = 0(Ax(")). 

Moreover, 

lim -fr.olUs = 0 and lim = {B,A) 

n—¥C!0 ’ n—^oc 

where jj • jji 2 denotes the L^-norm on R with respect to Lehesgue measure and Vr,o G (r G {6,s}) is 
non-negative, bounded and continuously differentiable. 

The first condition on the volume density functions is intuitive. The second condition will become clear 
later; it will be used to bound the impact of market orders and limit orders placed into the spread on the 
state of the book. 


1.1.2 Event types 

There are eight events - labeled A, ..., H - that change the state of the book. The events A, ..., D affect 
the buy side of the book: 

A := {market sell order} B := {buy limit order placed in the spread} 

C := {cancelation of buy volume} D := {buy limit order not placed in spread} 

The remaining four events affect the sell side of the book: 

E := {market buy order} F := {sell limit order placed in the spread} 

G := {cancelation of sell volume} H := {sell limit order not placed in the spread}. 


®This assumption, which may be generalized, considerably simplifies some of the analysis that follows. 
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We will describe the state dynamics of the n:th model by a stochastic process {S^'^}keK that takes values 
in the Hilbert space 

E :=M.xRx X L'^. 

(n) 

The first two components of the vector S^. stand for the best bid and ask price after k events; the third 
and fourth component refer to the buy and sell volume density functions relative to the best bid and ask 
price, respectively (visible and shadow book). We define a norm on E by 

||a||B := |ai| + |a2| + llaallis + ||a4||L2, a = (ai, 02,03,04) S (4) 

In the sequel we specify how different events change the state of the book and how order arrival times and 
sizes scale with the parameter n G N. 

1.1.3 Active orders 

Market orders and placements of limit orders in the spread lead to price changetQ. With a slight abuse of 
terminology we refer to these order types as active orders. For convenience, we assume that market orders 
match precisely against the standing volume at the best prices and that limit orders placed in the spread 
improve prices by one tick. The assumptions that market orders decrease (increase) the best bid (ask) price 
by one tick while limit orders placed in the spread decrease (increase) prices by the same amount have 
been made in the literature before and can be generalized without too much effort]^ However, this would 
unnecessarily complicate the analysis that follows. 

Thus, if the fc:th event is a sell market order (Event A), then the relative buy volume density shifts one price 
tick to the left (the liquidity that stood I tick into the book now stands I — 1 ticks into the book), the best 
bid price decreases by one tick and the relative sell volume density and the best ask price remain unchanged. 
Since the relative volume density functions are defined on the whole real line, the transition operators 

Ti"\v)i-) = vi- + Ax(”)), T^"\v)i-) = vi- - Ax^”^) 

are well defined and one has that 

= i-F {vti) (■), 41.(-) = 4:‘2(-) 

and 

4a=4">. 

The placement of orders into the spread will be modeled using the shadow book. If the /c:th event is a buy 
limit order placement in the spread (Event B), the relative buy volume density shifts one price tick to the 
right (the liquidity that stood one tick above the best bid in the shadow book now stands at the top of the 
visible book), the best bid price increases by one tick and the relative buy volume density and the best ask 
price remain unchanged: 

45+i(-) = (45) (■), 45+.(-) = 45(-) 

market order that does not lead to a price change can be viewed as a cancelation of standing volume at the best bid/ask 

price. 

market order whose size exceeds the standing volume at the top of the book and would hence move the price by more 
than one tick, is split by the exchange into a series of consecutively executed orders. The size of each such ‘suborder’, except 
the last, equals the liquidity at the current best price. Thus, by definition, a single market order cannot move the price by more 
than on tick. 
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and 


b<5. = b‘"> + a*("), 4”+’, = 4"’^ 


Remark 1.3. Notice that market order arrivals and limit order placements in the spread are “inverse 
operations”: a market sell order arrival followed by a limit buy order placement in the spread (or vice versa) 
leaves the book unchanged. 


1.1.4 Passive orders 


Limit order placements outside the spread and cancelations of standing volume do not change prices. With 
the same minor abuse of terminology as before, we refer to these order types as passive orders. We assume 
that cancelations of buy volume (Event C) occur for random proportions of the standing volume at random 
price levels while limit buy order placements outside the spread (Event D) occur for random volumes at 
random price levels. The submission and cancelation price levels are chosen relative to the best bid price. To 
guarantee convergence of the dynamics of volumes to the solution of a PDE, volumes placed and proportions 
canceled are scaled in a particular way. 


Assumption 1.4. For each k gN there exist random variables uj^,ujjf taking values in (0,1), respectively 
[0, M] for some M > 0 and random variables taking values in [—M,M] such that, if the k:th event 

is a limit buy order cancelation/placement, then it occurs at the price level (j G Z) for which 


IT 


C,D 

k 




The volume canceled, respectively, placed is 

oj^ • respectively wf* • 

Here value of the volume density functions at the cancelation price level and Au^"^ is a 

scaling parameter that describes the impact of an individual limit order arrival (cancelation) on the state of 
the order book. 


Some comments are in order. First, by construction, buy order placements and cancelations take place at ran¬ 
dom distances (positive or negative) from the best bid price. Since < M no orders are placed/canceled 

more than M/Aa;*-"^ ticks into the book. Volume changes take place in the visible or the shadow book, de¬ 
pending on the sign of of tt^. If >0, then the visible book changes; if tt^ < 0 the placement/cancelation 
takes place in the shadow book and the impact of the event on the state of the visible book (i.e. at price 
levels above the best bid) will be felt only after one or several price increases. 


Example 1.5. Suppose that the liquidity standing one tick above the best bid price after (k — 1) events is 
k-th event is an order placement of size Luj)Av^^'^ in the (buy-side) shadow book one tick 
above the best bid, i.e. tyJ) G [—Ax^^^O) and that the (k -\- l)st event is a buy limit order placement in the 
spread. Then Aa;*-"^ and the liquidity standing at the top of the book at (size of the 

order placed) is v^^f^ Luj) Av^^'^. This is how the buy order previously placed in the shadow now becomes 
part of the visible book. The role of the shadow book is further illustrated by Figure\M 


Second, the liquidity available for buying j GN ticks into the book after k events is given by by the integral of 
the volume density function over the interval [a;j"\ a:^"\), i.e. equals w[^^(a;^"^) • Ax^^\ When a cancelation 
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occurs at this price level, then the new volume is = (1—)''^b"fc (3:")-Aa:^"\ 

Hence cancelations are proportional to standing volumes. On the level of the volume density functions, 
Assumption 11.41 implies that 


4”fc+i(-) = (•) - (•), where {x) = ^ ^^(x). 

j^-OO 

Volume placements are additive. If the order a limit buy order, then the volume density function changes 
according to 

A (^) 

+ where := ^ )j(x). 

j=-oo 

In either case, the bid/ask price and standing sell side volume of the book remain unchanged: 

AVi(0 = »a('), Ba = d”'. 4i=4"’. 

Similar considerations apply to the sell side with respective random quantities and . 



Figure 3: Top left; Initial buy (red) and sell (blue) volume densities. The sell-side shadow book is uncolored. 
Top right: Second state after a sell-order has been placed in the shadow book at price level A{t) — Ax^’^\ 
Bottom left: Third state after a subsequent sell limit order has been placed in the spread. Bottom right: 
Third state if, instead of a subsequent sell limit order in the spread, a market buy order arrived. 

Assumption 1.6. For I £ {C,D,G,H} the sequences {w^lfegN and {7r|}fegN are independent sequences of 
i.i.d. random variables. Moreover, the random variables 7r| have C"^-densities with compact support. 

Lipschitz continuity of implies ||/(■±Aa;*■"^) —/(•)||oo = C>(A*^"^) as well as existence of a constant K < oo 
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such that 


ttI e [xf\xf^J\ -P[7r^ e [x^^\,xf'^)\ 


< 


<K 


\f{y + Ax^-^)-f{y)\dy 


/x(") 

J-1 


(5) 


Moreover, if we denote by the (scaled) expected changes in the heights of the density volume functions, 

due to cancelations and order placements, i.e. 


then 


.(1 . 

^ / f^{x)dx, 

3 

(6) 

-- 0(Aa;(")). 

(7) 


1.1.5 Event times 

The dynamics of event times is specified in terms of a sequence of inter-arrival times whose distributions 
may depend on prevailing best bid and ask prices. 

Assumption 1.7. Let {</Jfc}fegN be a sequence of non-negative random variables that are conditionally inde¬ 
pendent, given the current best bid and ask price: 


In the sequel we write (p{Ak, Bk) for ipk to indicate the dependence of the distribution of ipk on the best bid 
and ask price. Similar notation will be applied to other random variables whenever convenient. 

In the n:th model, we scale the time by a factor More precisely, we assume that the dynamics of the 

event times in the n:th model is given by: 


An) 

'fc-l-1 


= T, 


(n) 


where 4"') := 4“>) ■ At<- 


( 8 ) 


1.1.6 Event types 

Event types are described in terms of a sequence of random event indicator variables {(fk} taking values in 
the set {A,...,H}. We assume that the random variables 

<Zfe = <z(4”^4”^) (fcGNo) 

are conditionally independent, given the prevailing best bid and ask price and that their conditional proba¬ 
bilities 

pin) ^ pin)J^Bin)^^in)^ ^ 

satisfy the following condition. 


10 





Assumption 1.8. There are bounded continuous functions with bounded gradients —>■ [0,1] and a 

scaling parameter Ap^'^') —>• 0 such that 


•/(•,•) for / = A,B,E,F 

(1-Ap(")) ■/(•,•) for / = C,D,G,H 
1 
1 

Remark 1.9. The preceding assumption implies that an event is an active order with probability and 

a passive order with probability 1 — Ap^'^'i, independently of the state of the book. Conditioned on the order 
being active or passive, it is of type I with a probability p^{•, •) that depends on the current best bid and ask 
prices. We allow the above probabilities to be zero in order to account for the fact that no price improvements 
can take place when and to avoid depletion of the order 6ooA:@. 




pC+pD^pG_ 


-P^ = 

-p^ = 


The expected impact of each active order on the state of the book will be of order Aa;^”); that of a passive 
order of order Av^'^\ Because active orders arrive at a rate that is Ap*^"Atinies slower than that of passive 
orders, the relative average impact of active to passive orders on the state of the book will of the order 
^ scaling limit requires to equilibrate the impact of active and passive orders. In order to 

guarantee that there will be no fluctuations in the standing volumes in the limit as n —^ oo we also need a 
minimum relative frequency of passive order arrivals. This motivates the following assumption. 


Assumption 1.10. The scaling constants Ap^"^, Ax^^\ Au^"^ 


Ax(") • Ap(") . Au(" 


= Cl, 


for some a G (^, l) and constants Co,Ci > 00 


and At^^'i are such that: 
Ap(") 

and —. -^ = 0(1) 

(AtWj 


1.1.7 Active order times 

The previous two assumptions introduce two different time scales for order arrivals: a fast time scale for 
passive order arrivals, and a comparably slow time scale for active order arrivals. Inter-arrival times between 
passive orders are of the order At^") while inter-arrival times between active orders are of the order 
In order to see this, let us denote by ct). the arrival time of the fc:th active order. The number r\f.^ of 
events one needs to wait until the (fc -I- l)-st active order arrival can be viewed as the first success time in 
a series of Bernoulli experiments with success probability Ap(") and expected value The (fc -|- l)-st 

active order arrives at time 

(") _ («■) , An) . A (rj) 

^k+1 ~ ^ 

where 

^For simplicity we assumed that the the initial volume density functions vanish outside a compact price interval. Hence 
there is a positive probability of depletion unless one assumes that no further buy/sell side price improvements take place if 
the distance of the current best bid/ask price from the initial state exceeds some threshold 

^^For the results that follow, we will assume that '^■Ap^ ^ ^ ^ = 1 as n —)■ oo. Any other constant would 

require further constants in the limiting dynamics. 
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Since the random variables ip p m are conditionally independent and identically distributed, 

and are independent sequences. Because the conditional expected value 

m of 

given the prevailing bid and ask prices is independent of n S N. We assume that 
the mapping m(-, •) is Lipschitz continuous. 

Assumption 1.11. The conditional expected value m{B,A) of depends in a Lipschitz continuous man¬ 
ner on the prevailing pair of bid and ask prices {B,A). 


1.1.8 State dynamics 

We are now ready to describe the full dynamics of the state sequence. To this end, we put 

cl") _ / g(") An) An) (n)\ 

’ V- ’ ^b,k > ^s,k ) 

In terms of the indicator function 1^ := (^1 a {(t>k{B^k^, A ^^^,..., 1h {f)k{B^k^, A^^"' 

namics of the state sequence {<5'^”^} is of the form 

sh’i = sl”’ + p'”>(sl"’) 

in) 

if we dehne the random operator Vj,:E^Ehy 

where the matrix (5'^,"'^) equals 

— 0 0 


'))) the dy- 


( 9 ) 
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0 




0 0 

0 0 


V 0 


0 


0 


0 




’e,k Ax(") 




Here, the entries referring to shifts in the volume density functions, due to best bid and ask price changes, 
are given by 


jYfiAA rp{n) (An)\ _ (n) J\,An),E „(n) / _ A 

yb,k) ^b,k> ■“^+ \'^s,kj ^s,. 


M, 


W-" :=Ty> („<">)-„™, Mr' := T™ «I!) - »™ (W) 


-»S, := t 1"> («'”!) 




and the entries referring the volume changes, due to placement and cancelation of volume, are given by 

CX> 


J = -CXD 

Observing the dynamics in continuous time, we define 


S^-Ht):=si-^ and := for tGlrt\ril\). 


( 12 ) 
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Remark 1 . 12 . Overall, the state and time dynamics of our models are driven by the random sequences 
{(fik} (event times), {((k} (event types), {tt^} (placement/cancelation price levels) and ujI (placed/canceled 
orders). The joint dynamics of all models can be defined in terms of suitable independent families 


■= {{Tk{B,A),()k{B,A))(^A,B)m^: } {k € Nq) 

of independent random variables. In particular, the process {(5'*-"^^), S N) is adapted to 

the filtration 


4 ”) :=a 



: 0 < s < [ 



1.2 The main result 


Our main result is Theorem 11.131 It states that - with our choice of scaling - the order book dynamics can 
be described by a coupled ODE:PDE system when n —>■ oo: the dynamics of the best bid and ask prices will 
be given in terms of an ODE, while the relative buy and sell volume densities will be given by the respective 
unique classical solution of a first order linear hyperbolic PDE with variable coefficients. 

Theorem 1.13 (Law of Large Numbers for LOBs). Let be the sequence of continuous time 

processes defined in Dill) and suppose that Assumptions \1.2\ andhold. Then, for all T > 0 there 
exists a deterministic process s : [0,r] —E such that 

lim sup — s(t)||£; = 0 in probability. 

ra^oo jgjo,T] 


The process s is of the form s(t) = 


lit) 

v{-,t) 


where ^(t) = 


bit) 

a{t) 


is the vector of the best bid and 


ask price at time t € [0,T] and v{x,t) = 


Vbix,t) 


^ Vsix,t) ^ 

t S [0,T] relative to the best bid and ask price. In terms of the matrices 


denotes the vector of buy and sell volume densities at 


Ai-) := 


p^i-)-p^i-) 

0 


0 

P^i-)-p^i-) 


, B{-,x):= 


-P^i-)rix) 

0 


0 

-P^i-)f^ix) 


(13) 


the vector 


c(-,x) 


/ P^i-)f°ix) 
\ P^i-)f^ix) 


(14) 


and the function m(-, •) that specifies the expected waiting time between two consecutive active order arrivals, 
the function 7 is the unique solution to the 2-dimensional ODE system 


d'ljt) 

dt 

< 


7(0) 


^ Ah(t)) / 1 


Bo 

^0 


t€[0,T] 


(15) 
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and (vb^Vs) is the unique non-negative bounded classical solution of the PDE 


'^t{t,x) = m( 7 (t)) {l{t))vx{t,x) + B{-i{t),x)v{t,x) + c {t), x)"^ , (t,x) € [0, T] X R 

^(Oja;) =vo{0,x), a; € R 


( 16 ) 


The analysis of the limiting dynamics can be simplified by separating the randomness on the level of order 
arrival times from that of order types as shown in the following section. Subsequently, we give an explicit 
solution to the limiting PDE. 


1.2.1 State and time separation 

For the continuous-time process we write 

where G R^ describes the dynamics of bid and ask prices and Sv”^\t) G L^(R) x T^(R) describes the 
dynamics of the buy and sell volume density functions. According to the following proposition the process 
can be expressed as the composition of a state process and a time process The proof follows from 

straightforward modifications of arguments given in Anisimov [H p.l08]. 

Proposition 1.14 (State and time separation). The process 5^"^ can be expressed as the composition of a 
random state process 

and a random time process as 


5(”)(t) =??(”) (//(”)(t) - . 


The state and time process is given by 



r;(")(t) := 

for t G 

An) An) \ 

’ ^fc+1 J 

(17) 

k ■ and 

y(")(«) :=Tt^ 

for u G 

‘k ’ ‘k-\-lJ ’ 

(18) 

The time-change 

is then defined 

in terms of y^’^'i as 



p^^Ht) := inf{M > 0 : y^^\u) > t}. 


The advantage of the state and time separation is that the processes and can be analyzed separately. 
In fact, we will show convergence in probability 


lim sup 

n->-oo tg[Q T'] 


V^-Ht)-r,{t) 


= 0 
E 


and 


lim sup 

n-!-oO(g[g^y] 


p^-Ht)-pit) 


= 0 
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to limiting processes 77(t) = r/y(t)) and Since the state sequence takes values in the Hilbert space 

E, the time change theorem as proved, e.g. Billingsley [51 P- 151], then implies that 

lim sup — 77(77(t))||B = 0 in probability. 

n-)-oo tg[Q r] 

In our model, bid and ask prices are sufficient statistics for the evolution of the order book. In particular, the 
limiting behavior of the sequences and can be analyzed without reference to volumes. In Section [5] 
we prove the following proposition. 

Proposition 1.15. Let 7 be the unique solution to the ODE 


^ = Am)) 


7(0) = 


Bn 


, te(o,T] 


(19) 


Then 


lim sup 

rn-oo 0<t<T 




= 0 in probability. 


Moreover, the sequence of processes satisfies 


lim sup 

n^oo Q<:t<T 


p^'^\t) — p,{f) =0 in probability. where p, ^{f) = 




Once the limiting time-change process p has been identified, what remains to finish the proof of Theorem 
11.131 is to establish convergence of the volume processes 77),”^ to their deterministic limit. This will be 
achieved in Section |3| where we prove the following result. 


Proposition 1.16. Let u be the unique classical solution of the PDE 


Then 


ut{x, t) = A (j(t)) Ux(x, t) + B{x, j{t))u{x, t) + c (x, j(t )), (x,t) e R X [0, T] 

it(a;, 0) = vo(x), x G R 


lim sup 

n->oo 0<t<T 



u(t, •) 


= 0 
L2 


in probability. 


( 20 ) 


1.2.2 Explicit solution 

The PDE system (1501) is coupled only through the limiting price dynamics. In particular, the equations for 
the buy and sell side can be solved independently. For the buy side PDE, we can write 

=Ab{t)^+Bb{t,x)ut + Cb{t,x) 

[u{,(0,x) = Vb,oix) 
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where Ab{t) = Bb{t,x) := -p^(x), Cb{t,x) := p^{'y{t))f^{x). Using the 

method of characteristic curves, the PDE reduces to a family of ODE:s; see [H Chapter 3] for details. The 
characteristic equations for our buy side PDE read 


I ^ = -Abir) I ^ = Mr, x(T))ub + Cb{T, x{t)) 

la;(0)=C \ubi^,0 = '<Jb,oiO- 

The solution of this system of linear ODE:s as a function of the state ^ G K can be given in closed form: 


x{t,^) 


[ Ab{t)At 
^0 

exp ^ y Bb{u,x{u,^))duj(^Vbfi{0 + J exp ^ - y Bbiu,x{u,^))dujcb{s,x{s,^))ds^ . 


It describes the surface {(t,C) ^ Ub{t,x{t,^)) = Ub{t,^) given Ub{0,^) = fo.feCO}- The solution to the buy side 
PDE can be recovered from the solution to the ODE-system through 


Ub{t,y) = Ub{t,y + 

Due to our smoothness assumptions on the volume placement and cancelation functions it is not hard to 
verify that the solution is uniformly bounded with uniformly bounded first and second derivatives with 
respect to the time and space variable. Moreover, since the function Vbp vanishes outside a compact in¬ 
terval i Assumption 11.21) and no orders are placed or canceled beyond a distance M from the best bid/ask 
price fAssumption 11.41) . the function Ub{t,-) vanishes outside some compact interval I{T) for all t G [0,r]. 
Altogether, one has the following result. 



Proposition 1.17. Under the assumptions of Theorem \ 1.131 the PDE D61) has a unique solution Ub- The 
solution is uniformly bounded, with uniformly bounded first and second derivatives with respect to both vari¬ 
ables and there exists an interval I such that Ub(t, a;) = 0 for all t G [0, T] and x ^ I. 


1.3 Examples and applications 

In this section, we discus examples where the limiting dynamics can be given in closed form as well as possible 
applications of our order book model to portfolio liquidation problems. The situation is most transparent 
when orders arrive according to a Poisson dynamics. Let us therefore assume that in the benchmark model 
n = 1 orders arrive according to independent Poisson processes with smooth rate functions A"^(-, •),..., M{■, •)■ 
Then {<pk} is simply the minimum of the inter-arrival times of the Poisson processes. Standard arguments 
yield: 

To simplify the analysis we normalize the arrival rates so that rn{-, •) = 1. In the n:th model we scale the 
arrival rates of passive orders by and those of active orders by Ax^'^\ The arrival rates in the n:th 

model thus satisfy: 

•) = for / = A,B,E,F. 

and ^ 

A(")d(.,.) = 444 for / = C,D,G,H. 

^ ^ AtU) 
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Figure 4: Bid/ask price dynamics of Example 1 1.1 81 


Then the probability 


of the next event being an active order is of the order 




and: 




_ 

_A^M_ 

AC(-,-) + A^(-,-) + AG(.,.) + A^(-,-) 


for / = A,B,E,F 
for / = C,D,G,H. 


1.3.1 Price dynamics 

Within the Poisson framework we can readily set up the ODE for the limiting dynamics of the best bid and 
ask price process. Of course, the spread &(t) := a(t) — b(t) should be non-negative at all times. This can 
easily be achieved if we require p^(-) = P^{-) = 0 for 6(t) = 0. For most application it would in fact be 
sufficient to assume that the price dynamics depend on the prevailing best bid and ask prices only through 
the spread. If the stationarity condition 

A^ - A-® = A-® - A-^ = 0 

holds, then the spread is constant, &(t) = 6(0), and always positive if the we require the limiting initial 
spread to be positive. The following is a simple example where the spread is initially positive, then becomes 
zero, and eventually opens again. 

Example 1.18. Assume that the order arrival rates take the form 

( X^{b{t),a{t)) = A^(6(t),a(t)) 

\ X^{b{t),a{t)) - X^{b(t),a{t)) = 2y/a{t) - Bq, t € [0,T] 

Then the bid price is constant, b{t) = Bq, and a(t) = Bq + (t — for t £ [0,r]; see Figure 4- 

The probability of a price increase, respectively, decrease can also be approximated in terms of the arrival 
rates. If the next event is an active order arrival on the buy side, then 

hm = B-{t-) - hm = B^{t-) + 

n—fco A + A n—fco A + A 
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The limiting probabilities that the next price change is an increase/decrease of the best bid/ask price as well 
as unconditional probabilities of price changes can be computed analogously. Clearly, these probabilities 
increase in the respective rates. A further possible application of our model includes estimations of the 
expected time-to-fill of a limit order in the original model. This question has been studied by, e.g. Lo et al 
[25] . A simple approximation in our model framework is given by the first time the limiting price process 
hits the placement price level of the limit order. 


1.3.2 Volume dynamics 


We now turn to the volume dynamics. When the limiting price dynamics is constant, then the explicit 
solution to the bid-side volume density PDE given in (1221) simplifies to 


Ub{t,y) = ivbfiijj) 


Cbjy) 
Bb{y) L 




(23) 


where we write Bb{y) and Cf,(y) for Bb(^{Q),y) and Cb(7(0),y), respectively. From this we see that the 
stationary solution is 

(f \ _ f^(y) 

Bb{y) ~ 

In particular, the standing volume at the price level y increases in the order arrival rates and the “probabil¬ 
ities” /^(y) with which an order is placed at that particular level when a placement occurs. Analogously, 
they decrease in and /^(y). The following is a simple example where the spread is constant, but the 
prices are not. 

Example 1.19. Let us assume that 

X^{b, a) - X^{b, a) = X^{b, a) - X^{b, a) = 1. 

Then b(t) = bo+t and a(t) = b(t) -I- 6(0). Assume moreover that bo = 0 and that the buy side passive order 
arrival rates depend only on the spread: 

X^{b, a) = X^(b, a) = a — b. 

Thus, p'^’^{j{t)) = p^'^ ■ Let us further assume (ignoring the fact that our density functions need to be 
defined on compact intervals for simplicity) that 


f= ^ 


and / (x) = -fj-. 

pL, pU 


Then, 


Bb{l{t),y) =-I and Cb{'y{t),y) = e 
We compute from the general solution formula: 

Ub{t, x) = e“‘Mf,(0, x) + e“*'(l — e“‘). 
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1.3.3 Portfolio liquidation 


We close this section with a brief discussion of how our LOB model could be used to obtain endogenous 
shape functions for portfolio liquidation models. In models of optimal portfolio liquidation under market 
impact the goal is to find trading strategies that unwind (or purchase) a large number X > 0 of shares 
within a pre-specified time window [0,T] at minimal cost. It is typically assumed that prices are continuous 
(as in our limiting model), and that the expected distribution of the standing buy (or sell) side volume can 
be described in terms of a possibly time-dependent shape function 


/ : [0, T] X R R+ 


The benchmark case studied in the seminal paper of Almgren and Chriss [3] corresponds to a block-shaped 
order book where f{t,x) = S for some d > 0; discrete-time liquidation problems with more general shape 
functions have been studied in, e.g. [2]. For a given shape function a sell order of size Et submitted at time 
t € [0,T] moves the best bid price by an amount Dt defined through 

pDt 

Et= f{t,x)dx. 

Jo 

If we denote by •) the inverse of the anti-derivative of the shape function and if we assume for simplicity 

that there is no permanent price impact and that the cost of trading is benchmarked against the mid quote, 
then the cost c{t,Et) of trading Et shares is half the spread plus the market impact costs (see, e.g. [2] and 
references therein for details): 

c{t,Et)=Et + F-\t,Et)^ . (24) 

If orders can be submitted at discrete points in time (n = 1,..., N), the resulting optimization problem is 
given by: 

N 

min c{tn,Et^) s.t 


N 

n—l 


= X. 


(25) 


Our model can be used to derive shape functions from order arrival dynamics. Within the framework of 
Poisson arrivals, the buy side dynamics in the n:th model is equivalent in distribution to modeling passive 
order arrivals and cancelations at the price levels {i ■ Aa:^"^}igz in terms of independent Poisson processes 
with respective rates 


(n),C/ \ 

m (•) 

= A(’^)’^(-) 

' / 

) 

f^{x) dx 

{iez) 

m (•) 


p{i-\-l)-Ax^^ 

) / 

J i-Ax('^') 

) 

f^{x) dx 

(^ e z) 


Placement and cancelation rates at different price levels (relative to the best bid/ask price) and average 
volumes placed/cancelled can readily be estimated from flow data; the same applies to aggregated (over all 
price levels) arrival rates as well as market order arrival and spread placement rates and 

volumes in the spread. We refer to Massey et al. m for methods and techniques for estimating parameters 
of (non-homogeneous) Poisson processes and to Hautsch [501 HI] and references therein for applications of 
Poisson and general point processes to order dynamics and high-frequency trading. 
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Remark 1.20. A possible dijficulty would be to calibrate the volumes placed in the spread. In our model a 
spread placement is the result of many small events taking place in the shadow book while in practice this 
is one single event. On the other hand, one could use a low-dimensional parametric model for the initial 
state of the book (e.g. block-shaped) and the density functions f^{x) for I G {C,D} and x < 0 and try to 
calibrate the resulting spread placements to empirical data. Another possibility would be to work with the 
stationary solution; in this case the density functions are only required on the positive half line (and could 
hence be estimated from flow data). This is what is often done indirectly in portfolio liquidation models 
when one assumes that the (benchmark) price is a martingale and can hence be treated as a constant in the 
optimization problem. 


Once the placement, cancelation and active order arrival rates rif^\ ^^ and have been esti¬ 

mated (possibly as functions of the best bid/ask price), functions f^ and f^ that satisfy the above equations 
can be constructed, e.g. by applying a method outlined in Schmidt et al. [32] and the limiting model can 
be set up. 


Example 1.21. Let us assume that the tick size is and that we are given empirical arrival rates which 

satisfy 


(n),C ('n),D 


(n),C 

Vl+l 


{n),D 

Vl+1 


for all i G and some constants > 0. 


Then, the functions f^’^ are recovered as: 


Assuming that , the stationary solution is of the form Ub{t, x) = Kie for ki,K 2 > 0 


and 


Ub(t, x) dx 


— — K 
K2 


Let us now assume that functions f^’^ have been constructed along with the rates A"^,...,A^. Plugging 
the dynamics of the buy side volume density function into the portfolio liquidation problem, the resulting 
optimization problem reads 


N 


min 




-6(t„) -I- Uf, {tn,Et„) 


N 

s.t. ^ = A 

n—1 


where Uf, ^{t;-) denotes the inverse of the anti-derivative of the volume density function Ub{t;-). In the 
framework of Example 1 1.21 1 we arrive at the logarithmic impact function Ufl^{t, y) = ln(l — ^). 

Remark 1.22. If one assumes that continuous trading is possible and restricts the set of admissible trading 
strategies to absolutely continuous ones, then the liquidation problem can be stated in terms of the volumes 
Ub{t,0) {0 < t < T) at the top of the book. If both absolutely continuous and block trading is allowed as in 
\2^ . then the dynamics of the full book is again required. 


2 Convergence of bid/ask prices 

According to Proposition 11.141 the process can be represented in terms of a composition of a state 
process 77^”^ that jumps at deterministic times and a time process that accounts for the random 
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event arrival times. Prices change less frequently at times This suggests to introduce a second time 

scale - which will be referred to as active order time - defined by 

:= k ■ 

(n) 

along which to scale the price process. In order to make this more precise, let us denote by 2?^ the restriction 
of the operator to the price component of the state sequence and put 




V 


(n) 

■y.k 


= E 


V 


(n) 

1,1' 



7 ,cr 


(n) 

k 


where the second equality follows from the fact that prices do not change between times and — 1. 
Furthermore, we introduce the family of continuous time stochastic processes r]iy defined by 

(t) := for t € 


where 


(^(n) 

■■=tl + 

Pl,k+1 

^(n) 

( 

Vo 

.. 

- 1, 


(26) 


The quantity fj^l describes the state of the price process after the fc:th price change. The following lemma 
shows that the process r ]^'^, evolving on the level of event time, and the process , evolving on the level 
of active order time, are indistinguishable in the limit when n —>■ 00 . 

Lemma 2.1. For any T > 0 and e > 0, it holds that 


lim P 

n—^oo 



t/ 



sup 

E 

^(n) \ 

> e 

0<t<T 


k^O 



= 0. 


Proof. By construction, the two sums ^ ^ have the same expected value for 


any t £ [0, T]: 


^k—O 


7 ,A: 


Lt/Ad")j 

E >^^ 3 = E ^'^li- 

k=0 k=0 


As a result, it is enough to prove that 








Lt/At(")j 


lim P 

n—¥C!0 

sup 

0<t<T 

A.-0 

e 

>2 

= lim P 

n—¥oo 

sup 

0<t<T 

k—0 

- 1 

V 1 CM 

A 


= 0 . 


The random variables 


Aa;(") 


fc = 0,...,LT/Aa;(")j, n£N 
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and 

fc = 0,...,LT/At(")j, neN 

form triangular martingale difference arrays in the sense of Definition lA.ll with respect to the filtrations 
(n)}fceN and respectively. A direct computation shows that they are uniformly L^-bounded. 

'^k 

Thus, it follows from Theorem I A. 31 that for all a > 

lim P 

n—>-oo 

as well as 

lim P 

n—¥C!0 

Choosing a £ ( 5 ,1) and multiplying the inequalities in the above probabilities by and respec¬ 
tively, proves the assertion. □ 




Let 7 be the solution to the ODE (fT^ and consider the discretisation := 7 (s[”^). The next lemma 
shows that the sequence of expected price processes 7 ^"^ defined by 


where 


:= for t £ 

:=7i”^+E[^S(7r)^ 


7o 


(n) 


(27) 


0 

(n) 

0 


converges uniformly to 7 on compact time intervals. The proof is standard; we give it merely for completeness. 
Lemma 2.2. For any T > 0 

sup | 7 (fo(<)- 7 (<)| = C)(Aa:(fo). 
te[o.T] 

Proof. Active orders change prices by one tick. Moreover, by Assumption ll.Sl the conditional probabilities of 
an active order being a market buy/sell order or limit buy/sell order placement in the spread are independent 
of n £ N. Hence 


E 




= Aa;fo)-|A(7l.”^)(^ J 


(28) 


Thus, the sequence 7 *-"^ defines a special case of the classical Euler scheme for the ODE (fTKll and hence 
converges uniformly to its unique solution, see e.g. Hairer et al. |191 Theorem 7.3], with rate Acc^"^. □ 
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We are now ready to prove convergence in probability of the bid and ask prices. 


Proof of Proposition 11.151 

a) We first consider the convergence of the state process 77 ^”^ and claim that 

lim sup —>• 0 , in probability. 

In view of Lemma 12.11 we can write 

=7i">(o)+ i: 

k^O 


(29) 


= ^W(o)+ ^ J2 


fc =0 


fc =0 


up to some random additive constant that vanishes almost surely uniformly in t S [0,r] as n —>■ 00 . 
Adding and subtracting the sequence yields (again up to a vanishing additive constant): 




< 


7(")(t)-7(t) 

E -7(")(t) 


+ 


+ 


k=0 






For the first term, we deduce from Lemma \2 ^ that 


sup | 7 (")(t)- 7 (t)| =C>(At(”)). 

tG[0.T] 


(30) 


For the second term we use the Lipschitz continuity of the event probabilities {■, •) in order to establish 
the existence of a constant > 0 such that: 


\t/Ax^'^'^ \ 




Lt/Ax'^^J 




E ® 


_^W(7) 

= 

E ® 


-E 


/i-O 




/i-O 







^ E 


E 




A;=0 

[t/Ax^^^j 


-E 




'(n) 


%,k - Tk 

The third term corresponds to the noise-term of the price process. For each n S N, the sequence 

Vk ■■= - E [^3(57^"^)], fc = 0,..., Lr/Ax(-)J 
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is a martingale difference sequence. A direct computation shows that 

2 


supE|2/^|2<C-(Ax(")) . 

n t- V / 


Hence, the law of large numbers for triangular martingale difference arrays lTheorem lA.3l and Corollary 
IA.4[1 implies 


lim P 


sup \J2yk\>^ 

0<m<T/Ax<.”-) 


= 0 , 

just as in the proof of Lemma l2.II Thus, using Lemma |2.II again, we see that 


(t) - m = it) - m < ■ L, • ^ 




^(n) ^(n) 

<,fc - Tk 


o(l) in probability 


for some additive term of order o(l) uniform in t S [0, T], As a result, (l29l) follows from an application 
of Gronwall’s lemma along with Lemma l2.II 

b) Let us now consider the cumulative “active order time process” 
y(")(t)= ^ C^") 


k^O 




= • 

1 . . ... 

= A.<"| y: ™ ( 5 !,^) + A.<"| . y: (c‘“> - E [c'”' 


\t / \ 




fc=0 


fc =0 


(31) 


By the above established uniform convergence of to 7 in probability and because the function m 
is Lipschitz continuous, the first sum converges to the function 


y{t) = / m{j{u))du. 
Jo 


(32) 


Applying the same arguments as above to the martingale difference sequences 

t 


c<"> (5 S)-e[c<"> (5S); 

we see that the second term vanishes uniformly in t G [0j2”] in probability. Thus, 




lim sup - y{t)\ = 0 in probability. 


Since and y are increasing functions, their inverses and y exist. By continuity 

sup in probability as n ^ 00 

te[o.T] 


and 


^k'it) = {y 


1 

y'iy~^it)) 


1 _ 1 


(33) 
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Since both the state and the time process converge, we conclude from the time change theorem that 


sup 

tG[0,T] 




—> 0 in probability n —>■ oo, 


where -j{t) = and 




AmKm / 1 \ Aj^jt)) / 1 

w(7(/i(t))) \ 1 / \ 1 


□ 


3 Convergence of volume densities 

In this section we prove Proposition [Td 6 l To this end, we denote by ’) ^^e restriction of the operator 


to L'^ X i.e the restriction of V]!'’ to the volume components of the state process. We need to show 
that the sequence {r]i^^}neN of x L^-valued step-functions defined recursively by 




Tyf") (t, •) := for t € [4”^, 44i) 


where 


(ra) 

%,k+l 


(n) 


(n) 

% 


_ „(") 


-v. 


(ra) 




(34) 


(35) 


converges in probability in to the unique solution of the PDE (1201) . We will show convergence in several 
steps. In a first step we find a convergent discretization scheme of the PDE which is coherent with the order 
book dynamics. Subsequently, we link this scheme to the expected dynamics of the volume densities. 


3.1 A numerical scheme for the limiting PDE 

For any n G N, the scaling parameters and define a grid {(4"^2;^"4} on [0,T] x R through 

= k ■ Atfo) {k G No) and = j ■ (j G Z). In a first step, we approximate the unique solution 

u : [0, T] X R — >■ R^ to (I20p through a sequence of grid-point functions : [0,T] x R — >• R^. To this end, 
we put 

' p^m)) 

0 


A (f) [ P^(7W) 0 


Ahit) := 


0 


p (tW) 


and 


E(t, x) := B( 7 (t), x), g{t, x) := c( 7 (<), x). 
Furthermore, we introduce operators 7^4^ that act on u G according to 

v!f'\v) ■.= V + Apfo^ • Aii(t) v(- + Aajfo^ ~ 

+ Apfo) • AL(t) k- - Axfo4 - z;(-)' 

+ Au(fo-(l-A4fo).[E(t,.)-u(-)+ff(t,-)]- 
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( 36 ) 


The sequence of grid-point approximations is then defined recursively by 

•) := for t e [4", 44i) 

where 


^k+1 


yin) 


=«:s,(d"’) 


(37) 


■‘0 “ ^0 

The sequence of step-functions essentially describes a discretized limiting volume dynamics of the 

order book. We benchmark this dynamics against the expected pre-limit volume dynamics when prices 
are replaced by the limiting dynamics. More precisely, we introduce another sequence of step functions 
: [0, T] X R M by 

•) := 4”^ for f€[4",4^i) (38) 


(") 


where 


“S. =4“’+e 


, for (£K, ()+’,) 

‘ ‘“^7(4"’),4"’) 




(39) 


(n) (n) 

K ■=yy- 


In a first step we are now going to show that the grid-point functions approximate the solution u of our 
PDE. Subsequently we show that the PDE can as well be approximated by the functions 

Proposition 3.1 (Convergence of the numerical scheme). Assume that the assumptions of Theorem 
hold. Then, the processes define a convergent finite difference scheme of the PDE i2(M . i. e. 

sup \\u^'^\t,-)—u{t,-)\\L 2 ^ 0 as n ^ oo. (40) 

te[o.T] 

Proof. For t G [0,T] the local truncation error associated with the grid-point approximations is defined 
as 

(41) 


C'^'^\t,x) := (u{t +At^'^\x)-'H[’^\u{t,-)){x)'j . 

Smoothness of the solution u (bounded with uniformly bounded first and second derivatives) along with 
Assumption 1 1.1 Ol implies that the following estimate holds uniformly in t G [0,r] and x G M: 

d''^\t,x) = + At^’^\x) — u{t,x) — • An(t){u{t,x + Ax^^^) — u(t,x)) 

— Ap^^l ■ AL(t)(u(t, X — ^)) 

— Ax^”^ • (1 — Ap^”^) • [F{t, x)u{t, x) — git, x)] ^ 

= (ut{t,x)At^'^^ -4 o(At("'4 - ^4"^ • Afl(Ma;(t,x)Ax("') -ko(Ax("'4) 

- Ap(") • Az,(-ii,,(t,x))Ax(”) -k o(Ax("))) 

— Au^"^ • (1 — Ap^"^) • [F{t, x)ult, x) — git, x)] ^ 

= ut{x, t) - {AL{t) - Anf)) U 3 ,(x, t) - F{x, t)u{x, t) - g{x, t) + o(l) 

= Ut{t, x) — A(t)ux{t, x) — Bit, x)u{t, x) — c{t, x) + o(l) 

= 0(1) (42) 
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The solution u{t;-) the PDE vanishes outside a compact interval for all t € [0,T] (Proposition [TTT7|) and 
hence so does •). As a result, we also have 

lim sup = 0. (43) 

n^oo 


From (j41|) . one has that u(t + •) = •)) + ■). In terms of the error-function 


this yields 


(t, x) := {t, x) — uit, x) (44) 

5S(")(tW ;.) = 52(-)(4;\;.) - ir(4"^ + •) 

= •)) - •)) - At(")/:W(4-\.) 

due to the linearity of the operator Using this property iteratively, and putting '■= o • • • o 


one finds 


<5m(44i> •) = (H(”4'=(5u(" 40, •)) - At(") ^(H(”^)'=-*(A(”4tr ). •)■ 


(45) 


2=0 


in) 

From the definition of 1-L\ together with the fact that the functions F and g are uniformly bounded by 
assumption, one finds (for large enough n G N): 

||h44^)IU^ < \\v\\l^ (l - ^p‘'^\AR{t) + As{t))) + ||Tf4u)|U2 . • An{t)) 

+ \\t[^\v)U. • (Ap(-) . AUt)) + Au(") (IlFlUlklU^ + ||5|U). 

Thus, by the isometry property ||T£”4 'c)||_l 2 = ||u||i 2 of the translation operators, there exists a constant 
C > 0 that is independent of t G [0, T] such that 


sup \\nl"\v)\\L 2 <l + CAt^^\ 
lhlL2=i 


In particular, since k < J 


sup 

k 


•)) < (l + CAt(’^))LT/AtWj||^^(n)(Q^ 


Il2 


< e^^||5u(")(0,-)IU= 

= o(l), 

where the last equality follows from Assumption 11.21 Similarly, from 

(^(n))fe-.(^(n)(^W)^.) =0(1). 


sup 

k,i,n 


L2 


Using the same arguments as before, one concludes that the sum in (1451) vanishes uniformly in time. Hence, 


lim sup —u{t,-)\\j ^2 = 0. 

n-)-oo tg[Q r] 


□ 
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Next, we show that the functions also approximate the PDE. More precisely, the following holds. 
Proposition 3.2. Under the assumptions of Theorem \l.lS\ 

lim sup •) — •)||i 2 = 0. 

rt-J-oo 

Proof. The proof is similar to that of Proposition 13.11 By analogy to the operator we introduce an 
operator on by 

n["\v) :=v + • Anit) fn(- + - n(-)' 

+ Ap(”) • AL{t) |'n(- - Ax(”)) - n(-) 

+ An(") • (1 - Ap(")) •) • v{-) + g^^\t, •) 


where for t € [4”\4+\): 






\t,-) ■■= 


/W'^(.)-p^(7(4”4) 


For the error function {■) := ~ then obtain 

= ni'^\Svi-\-)) + Au(") ■ (1 - Ap(")) . (<5 fW(.) . ui"4-) + 

=: ni"\Sui"\-)) + Au(”) • (1 - Ap(")) • £(4”V) 

where 

' -(/(").^(.)-/^(-))-p''(7(4”4) 0 


Sg^^Htr) ■■= 


0 J’ 

(/(n),D(.)-yD(.)) .p^(9(4-))) \ 

(/("^’^(•)-/^(-))-p^(7(4”4) r 


By construction, the grid-point functions are uniformly bounded. As it result, it follows from Assumption 
If.41 that 

lim sup ||£( 4 "'V)I|l 2 = 0. 

k=0,...,[T/AtM\ 

One can now proceed as in the proof of ProDOsition l3.ll to conclude that 

lim sup •) — •)||i 2 = 0. 


□ 
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3.2 Expected volume dynamics and discretized PDEs 


(n) 

To show the convergence of volume density functions we compare the random states ryi ^ with the determin¬ 
istic approximations of the limiting PDE obtained in the previous subsection. For this, we introduce the 
deterministic step function valued processes : 


where 


•) := 4”^ for t e [4”>4+i) 






V' 


,k\Pl,k^^k ) 


t(") 




(46) 


(47) 


The process describes the expected dynamics of the volume density functions for the actual price 
process; in particular, is a stochastic process. By contrast, the process describes the dynamics of 
the expected volume density functions when the random evolution of bid and ask prices is replaced by its 
deterministic limit. We have: 

.) - u{t, .)||^. < ||^(")(i,.) - 5l(-)(t, .)||^. + •) - u(")(t, ■)U. 

+ ||w(-)(t, •) - Oil + Wu^^Ht ,.) - u{t, .)|U2. 

The last two terms are deterministic and converges uniformly to zero by Propositions 13.1] and l3.2] It remains 
to show convergence of the first two (random) terms. This will be achieved in the following two subsections. 


3.2.1 Estimating the price impact of expected volume dynamics 

The term •) — OIIl^ measures the impact of the noise in the price process on the expected 

standing volume. The following proposition shows that it converges to zero almost surely (i.e. for almost all 
price processes), uniformly over compact time intervals. 

Proposition 3.3. Under the assumptions of Theorem \l.lS\ it holds that 

sup ||u^"^(t, •) — OIIl^ —> 0 in probability as n ^ oo. 

te[o,T] 


Proof. We argue again as in the proof of Proposition 13.II Analogously to the operator defined in the 
proof of Proposition 13.21 we define for t S the operator 

:= ?; -I- • A^^\t) v(- + — v(-) 


where for t € [4"\44i): 


+ • AL(t) k- - Ax(”4 - v(-) 

+ -(l-Ap^")) [F(”)(t,-) •u(-)+4”4i,-) 


:= 


0 
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and •) and •) are defined analogously. Let us further put := 

for I = A, B, E, F and 


SA 


(n) 

R,k 


6pl^^’^ 0 

0 Spi^^’^ 


SA^lI ■■= 


Spi^^'^ 0 \ 

0 Spi^^’^ ) 


and denote by <5/^,"^ and 
Then the error function 6u^^\-) 
the operator as follows: 


the corresponding quantities for cancelations and limit order placements. 
■= ~ satisfies = 0 and can be represented in terms of 




(n) 




where 


4”'(4”>,.) := 


i,(„) 


,(") 




+ 


Ap(") 

At(") 


<5A 


(n) 

L,k 


,(") 






Corollary IB. 41 establishes 


< L and < L ■ Ax^"^ 

for some constant L < oo that is independent of {n,k). Using our assumptions on the placements and 
cancelations functions and event probability functions along with the fact that 

lim sup “ 7(4"4| 0 in probability, 


this implies 

lim sup ||£[,”44”\-)IU= = 0- 

fe=0,...,LT/At(")J 

We can now argue as in the proof of Proposition 13.II to conclude. 


□ 


3.2.2 Convergence of volumes to their expected values 

In this subsection we apply a law of large number for Hilbert space-valued triangular martingale difference 
arrays (TMDAs) in order to establish the missing convergence to zero of the distance between 4"^ and . 
More precisely, our goal is to prove the following result. 

Proposition 3.4. Suppose the assumption of Theorem \l.lA hold. Then, 

sup \\g^^\t, •) — u^'^\t, OIIl^ —^0 in probability as n ^ oo. 
ie[ 0 .T] 


Proof. Using the definition of 44 ia (1551) and ia (H51) we see that 


t(-) = e4”^ 


v,k ’ 
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conditioned on the price process. As a result, 


Lt/Ad")j 




fe=0 


In order to establish convergence of the sum to zero uniformly in time we introduce the L^-valued triangular 
mart ingale- difference- array 

n” := -Diii; •) - (48) 

If we can show that there exists ^ such that 


/ 1 N 2/3 

then Theorem IA.3I and Corollary IA.4I would guarantee that 

m 


lim P 

n—^oo 


sup 


= 0 


0<m<lT/AtM\ 

and the proposition would be proved. To establish (HHl) . we need to bound the following terms: 


supE 

n,k 


1 p 

Dk,Hk v,k 


'-DkMk 


Av^ri) 

Ax(") 




supE 

n,k 


^Ck,Gk 


Aa;(n) 


M. 


(n),C,G (ri) 




Au(”) 


1 ““ (n) 

ICfc.Gfc Afy j, rj^ f. 


supE 

n,k 

supE 

n,k 





- "l 

^Vb,k J 

-E 



{^Sk) 

- ^ 

^Vb,k J 




- z?^”^ ^ 
,k J 

-E 

^Ek,Fk 



_„("■) V 

^Vs,k J 


• 

2 ■ 


L2 

1 

2 ■ 

. 

L2 

1 

2 

- 

L2 


(49) 


(50) 


This is done in Lemma lB.Gl in the appendix. In particular, this lemma shows that the placement/cancellation 
and shift terms of of the order 

C> ((Au(’^))2) and O ((Az;^"))^^-h (Az;^"))^-^) , 

respectively. Since a G (1/2,1) the assertion follows for /3 := minja, 1 — a/2}. 

□ 


4 Conclusion 

In this work a law of large numbers for limit order books was established. Starting from order arrival and 
cancelation rates for all price levels, we showed that the LOB dynamics can be described by a coupled 
PDE:ODE system when tick and order sizes tend to zero while arrival rates tend to infinity in a particular 
way. A key insight is that the scaling limit requires two time scales: a fast time scale for passive order arrivals 
and a comparably slow time scale for active order arrivals. The proof of convergence of volume densities was 
carried out in three steps: We first showed that the expected LOB dynamics resembles a numerical scheme 
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of hyperbolic PDEs plus noise, provided the random price dynamics is replaced by its deterministic limit. 
Subsequently, we showed that the impact of the noise in the price process on the volume dynamics vanishes 
in the limit. Finally, we used a law of large numbers for triangular martingale difference arrays to prove 
that the LOB model converges to its expected value. Our model allows for approximation of key order book 
statistics such as expected price increments, expected standing volumes at future times and expected time 
to fill. 

Several questions remain open. First, it would be interesting to establish a CLT or, more generally, a 
diffusion approximation for LOBs. Based on the idea of having different time scales for active and passive 
order arrivals, Bayer, Horst and Qiu have recently established a first SPDE scaling limit for order books. 
However, they assume that cancellations are subject to additive (rather than multiplicative) noise so volumes 
may become negative. Second, it would be interesting to solve models of optimal portfolio liquidation based 
on our limiting model. For this end, the model should ideally be calibrated to market data. 


A A Law of Large Numbers for Banach-Space-Valued TMDAs 

Here we prove a law of large numbers for triangular martingale difference arrays taking values in real separable 
p-uniformly smooth Banach spaces. A Banach space E is called p-uniformly smooth, where p G (1, 2], if 

Pe{t) = supjll^^ll + ll^y^ll - 1 : ||a;|| = l,\\y\\ = t| = ©(r^). 

All Hilbert spaces are 2-uniformly smooth by the parallelogram identity. The spaces C, and are not 
uniformly smooth. 

Definition A.l. A family of random variables k = 1,..., n, n = 1, 2,... defined on some probability space 
(H, P) is called a triangular martingale difference array (TMDA) with respect to a family {J^"}n=i, 2 ,... 
of filtrations, = 1)2,... the sequence y",...,y^ is an -martingale difference 

sequence (MDS), i.e. 

E [2/fe Vfe-i] = 0. 


If {yj^} is a TMDA, then for all n = 1, 2,... one has 

k 

j=i 

that is, partial sums are martingale. For such martingales, Pisier |30j proved the following moment estimate. 

Lemma A.2. Let E be a real separable p-uniformly smooth Banach space (1 < p < 2). Then, for all r >1 
there exists a constant C > 0 such that for all martingales 


with values in E, we have 



n>l 


E 


sup \Xn 

n>l 


< CE ^ 


r/p 


\n—l 
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This lemma allows us to prove the following law of large numbers for TMDAs. 

Theorem A.3. Let y^, k = n = 1,2,... he a TMDA taking values in a real separable p-uniformly 

smooth Banach space E for 1 < p < 2 such that 


supE|y^|P < oo. 

n,k 


Then, for all P > 0 such that (3 ■ p > 1 one has for all e > 0 that 


lim P 

n—^oo 


sup |^ 2 /fe|>e-n^ 


l<m<n 


= 0 . 


Proof. By Markov’s inequality 

m 

sup \ J2y^\>e-nP 

<m<n 

for all g > 1. Thus, it follows from Lemma [A .2 1 that 


< 




E 


sup I^Pfcl 


p 

m 

sup |^yfc|>e-n^ 


' n 


l<m<n 


.k^l 


-\ g/p 


< Cn 


-P-g+i 


for a generic constant C > 0 since the random variables y'f. have a uniformly bounded p:th moment. Hence, 
the assertion follows as soon as—/3-g+^<0, which holds for all g > 0 as /3 • p > 1. □ 


As an immediate corollary from the preceding theorem one obtains the following law of large numbers for 
TMDAs. 


Corollary A.4. Let k = l,...,n, 

smooth Banach space E such that 


n = 1,2,... be a TMDA taking values in a real separable 2-uniformly 
sup(n2%|y^|2) <oo 

n,k ' 


for some fi > \. Then 



in probability. 


Proof. We apply Theorem lA.31 to the TMDA 


''Ti Bn 

Vk ■= n^Vk- 


Then 


lim P 

n—^co 


sup 

l<m<n 


k=l 


= 0 . 


Hence the assertion follows from: 


m 


m 

sup 

= p 

sup \^yl\>c 

l<m<n 


l<m<n 


□ 
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B Properties of volume density functions 


In this appendix, we prove some properties of the volume density functions. In particular, we show that the 
sequences take values in almost surely. We first use an induction argument to establish a useful 

representation of the volume density function. 

Lemma B.l. The buy side volume density function satisfies: 



Proof. For k = 0, the equation holds by definition. Let us therefore assume it holds for all k < p. For 
k = p + 1 one then obtains 



□ 


B.l Boundedness of volnme densities 

Using the isometry property of the translation operator we deduce that the norm of the volume density 
function can be estimated from above by considering a model with only passive order placements. In a 
similar way we can estimate the expected order book hight at any given price tick. More precisely, we have 
the following result. 


Lemma B.2. The expected L'^-norm of the volume density function is uniformly bounded: 


sup < C 

■iGN, fe=0,...,LT/At(")J 


(51) 


for some constant C < oo. Likewise, the expected order book hight is uniformly bounded, i.e. if we put 
then : 

sup 

jGZ. nGN. k=0,...,lT/AtM\ 


7 L. I <C. 

'Vh 5 k I 


(52) 


Proof. It is enough to consider a model with only order placements where = 1 a.s. W.l.o.g. we may also 
assume that < 1 a.s. and = 0. Furthermore we may as well use a representation of the volume 
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densities in absolute rather than relative coordinates. In such a model, E||77^"^^||^2 is of the form 


^ ^ JGZ \i=l / 


where := 1 r m m and the distribution of wf properly adjusted to account for the representation 

in absolute coordinates. Since the random variables have compact support, only finitely many summands 
are non-zero and we may rearrange terms. Using conditional independence of the placement variables though 
time, this yields: 


nh 


in) 


2 fe 


Vb,k\\L^ 




i=l ji 


Ax(") 






Using the fact that no placements take place at price levels with a distance of more than 1 from the prevailing 
best bid/ask price: 


E 




I<1} 


Eal-^Ea^rJl < Wf^Wl^Ax^^^y 




In particular, the inner sums extend over at most + I terms. As a result our scaling assumptions 

guarantee that 

sup E||77i”^J|2 2 < oo. 

n,fc=l.[T/At(")J 

The second assertion follows analogously as 




Vi=l 


□ 


B.2 Norm estimates 

The next result will be used to prove a Lipschitz continuity property of the grid-point approximation of the 
limiting PDE. 

Lemma B.3. There exists a constant C < 0 such that for all n gN and k = 0 ,..., [T / At^'^'i \: 


E 


ri’-'irA) - lA] , <C-Ax("i. 


L2 
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Proof. Using Lemma FB . 1 1 and the linearity of the translation operator it follows that a.s. 




= T 


U) 




(n) 

-^v,,k 

E fc ~ 1 -I A 
i=0 




Eto If 


2 = 0 ^ 


. V ((r 

Aa;(”) I V + 


2 = 0 


E fc— 1 A 
j = t + l j 


4" i-S) - »S) 


E fc — 1 B 
.7 = 2 + 1 ■*■.7 


E fc — 1 1 ^ ^ 

.7 = 2 + 1 


T 


(n) 




-M, 


{n),D 


>f) 


(53) 

(54) 


) ([rW (m;”>-C,W) _ M<"54',W] If) (55) 


Taking the expected value and norms in (I55|) we find: 


E 


< 


-i(n) 


(iS) 

(*’) 


r(") /„(") A _ „(") 

\^vt,k) ^vt,k 

+) 
f),0 


Z,2 


Au(") 
L2 All"') 

Ar)(") 

Aa;W 


fc-i 


+ 


^ ||t|”^ (e ) - E 


i=0 

fc-1 


M. 


{n),D 


L2 




2 = 0 


^(n),C (n)' 

V,l 


-E 


fy^(»).C (n) 

U,2 >Vh,l 


L2 


By Assumptions II.21 and II.61 there exists a constant K < oo such that 


T 


(n) 


(+>) 


(rt) 


L2 


(56) 
< A:Ax(’") and 


AxU 


(e ) - E 






L2 




D 


L2 


< KAx^’^l 


As for the cancellation terms, independence of the event dynamics from the standing volumes yield: 


(e ) = r|”) (/(")’^Er?5), 


1 


Ax(") 


E 


(")] ^ f(")+E7?^"^- 

2 /Uh.2 IVh.l’ 


In view of the second assertion of Lemma lB.2l and using the fact that /(")+ is bounded along with Assumption 
[m we find a constant K < oo such that: 


Aa;(’^) 


T' 


(n) 


E 


M, 


(n),C (n) 


% 


-E 




(n),C (n)' 


L2 


< 

+ 

< K 


■> (/(">■=) e[tI"1(,w)-,w 

(e||t<”> (,(;■> )-,t>|U. + Ai<”>), 


Altogether, we arrive at the following estimate: 


E 


rpH /„U) A 

[^Vb.kJ 


(n) 

^v,.k 


[T/Ad")j 


L2 


< A:Ax(”) + KAv^"^ 


i=0 


E 


T 


(n) 

+ 




(n) 

^v,,^ 


L2 


(57) 
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Hence, it follows from Gronwall’s lemma that 


sup 

fc=0,...,T/At(") 


E 


T 


(n) 


yivb,k j 


(n) 


L2 


= o 




□ 


Corollary B.4. There exists a constant C > 0 such that 

■)\\l2 < C ■ 

Moreover, 

sup OIIl^ < oo. 

nGN,tg[0,T] 

Proof. In order to establish the first assertion we represent the functions as 


where 


= Ecl”' 

(58) 

'cffi -c1"’ + »S(7(4”').c1”>) 

(59) 

g(") j'_ 

•— <^0 



For t &[k ■ At("\ (fc + 1) ■ A*^")) the preceding lemma then implies: 




E 


-i(n) 


(cf) - c 


(n) 

k 


L2 


< CAa:(" 


The second assertion follows from (1551) together with Lemma IB. 21 

= iiEci”^iiL^ < E|ici”^iu. < c. 


□ 


Using Lipschitz continuity of along with the point wise shift estimate of the initial volume densities of 
Assumption 1 1.2 1 the following result can be established by analogy to Corollary IB.31 


Corollary B.5. 


sup 

iiGN, jeZ, k=0,...,lT/AtM] 


E 


(n),j+l _ (n),j 

'Vb,k '>Vb,k 


= 0{Ax^^^). 


(60) 


We close this appendix with norm estimates which are key to the proof of Proposition [XT] 
Lemma B.6. The following norm estimates hold: 


supE 

n,k 


supE 

n,k 


ICfc.Gfc 


1 

Dk,Hk ^^(n) 'I’.k 


-E 


Av^^‘ 




{n),D,H 




Aa;(") ” 




1 — iy^{n),C,G (n) 

ICfc.Gfc ^^(n) ^v,k 'tv.k 


supE 

n,k 

supE 

n,k 




(lA) 

- 'l 

^Vb,k j 

-E 




- "l 

^vt,k J 



(iS) 

- ^ 
,k J 

-E 

ISfc.Ffc 



- ^ 

^Va,k J 


= o 


= o 


O ((Au(”)) 2“ + (Au^”)) 


2-ot 


2-ot 
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Proof. The first two estimates follow from boundedness of the density functions along with independence 
of the event dynamics from volumes and (1511) . In order to establish the third and fourth estimate we have 
to prove that 


sup 

n,k 




'lvb,k\ 





To this end we use a representation of — ri^\ as in Lemma IB.3I but in absolute rather than 

relative coordinates. This means that the shift terms drop out of the representation but the probabilities of 
placements and cancellations need to be properly adjusted. 

Assumption 11.21 allows us to bound the impact of the initial condition (l53l) by a term of the order 
To compute the norm of the sum in (I54|) we need to compute a term of the form 





where 


<■)"’ = • 

f 1 

if Trf 

e [(j - 1) ■ Ax(”r j 


-1 

if Trf 

G [j ■ Ax(”r 0 + 1) 

Ax(”)) 


i 0 

else 



constant 

K < oo such that 


E| 

“Sir 

< 




< 


il<i} 


where the second inequality follows from ([5|); the indicator functions account for the representation of volumes 
in absolute coordinates. Using the fact that the random variables have compact support and that events 
are conditionally independent through time, we can now argue as in the proof of Lemma IB. 2 1 to deduce that: 




Aa;(”) 

Aa:(") 


jGZ \i=l / 


(Au(”))2 


= K 


Aa;(") 
((Aa;("))^ 


\2o! 


(Aa;("))2-“ 


(Aa:(’"))4 


To compute the norm of the sum in (1551) we need to compute a similar term, but with replaced by 




•St 

if Trf 



[ 0 

else 

(“Sir 

< KAx^^h 


from volumes along with (1601) and Lipschitz continuity of we also obtain a constant K < oo such that: 


Ea, 


(n) 


< 


< K 


e [j ■ (j + 1) • Ax("))] • Er;^’^' - PK^ G [O' - 1) . Aa:("), j • Aa:^"))] • 

(Aa:("))E G [0 - 1) ■ ■ Ax^"))]! . 
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Hence it follows from Corollary [B3] that 


Ea, 


(") 

id 


< K 




and the assertion follows as in the case of placements. 


□ 
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